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I. INTRODUCTION 

A unitary operator U is defined by W = mean- 
ing that its adjoint (Hcrmitian conjugate) is equal to its 
inverse, implying that U'U = UlP = I, where / is the 
identity operator. For discrete systems, unitary oper- 
ators can be represented as matrices. In general, the 
determinant of a unitary matrix is a unit complex num- 
ber, and thus they are classified as unimodular, mean- 
ing |det([/)| = 1, since the modulus of the determi- 
nant is of unit length. An important sub-group of the 
n-dimensional unitary group U(n) is the special unitary 
group SU(n), where "special" means det(U) = 1. 

Unitary operators have many uses in quantum mechan- 
ics. First and foremost, they describe all transformations 
that simultaneously preserve state normalization, Her- 
miticity, and non-negativity, thus preserving the physi- 
cality of quantum states, both pure and mixed. Time 
evolution of closed-system quantum states is a unitary 
process [1] . In quantum computation, unitary operations 
describe quantum logic gates [2, 3]. In particle physics, 
subatomic particles can be represented using special uni- 
tary groups [4-7] . Originally, unitary matrices were stud- 
ied to provide alternative descriptions for the compli- 
cated rotational motions of gyroscopes, giving rise to the 
well-known Cayley-Klein (CK) parameterization of rota- 
tions using complex numbers [1, 4]. 

A noteworthy feature of unitary operators is that any 
product of unitary operators is also unitary. A related 
and useful property is that any n-dimensional unitary 
matrix can be factored into a product of single-qubit uni- 
tary matrices [3], which provides a means for realizing 
any unitary operation experimentally, as shown in [8]. 
Due to their importance, there is a great wealth of infor- 
mation about unitary matrices, so in the present work we 
limit ourselves to only the bare essentials for our purpose. 

The goal of this paper is to present a canonical form 
for unitary matrices in any dimension, subject to the con- 
straints that it be as simple as possible without sacrificing 
generality, and that it be readily useful for common ap- 
plications. First, we give a general prescription in terms 
of CK parameters and constraints, then we incorporate 
the constraints through hyperspherical parameterization. 
Finally, we examine several physical applications. 



II. CAYLEY-KLEIN PARAMETERIZATION OF 
SPECIAL UNITARY MATRICES 

Before proceeding, we must establish some practical 
conventions. First, let our canonical form for U have 
only single terms in its first row and first column, where 
the top row and the left column are both first. The mo- 
tivations for this are as follows. As a change of basis, U 
transforms the first basis element |1) to a vector |1') with 
coefficients equal to the first row of U, such as 

[U hl t/i, 2 V|l)\ _ {U 1A \l) + U h2 \2)\ = (\1')\ 
\ U 2 ,i C/2.2 J\ \2)J ~ \ U 2A \l) + U 2 , 2 \2) J-\ |2') )■ 

(1) 

Additionally, if a pure state in the first basis state = 
|1) is operated on by U, it is transformed into a state 
whose coefficients arc the first column of U, such as 




The origin of this first convention is two-fold; in general 
it is only possible to get at most one row and one column 
of a unitary matrix to have single terms in all entries, 
and focusing on the first row and first column vectors 
(leading vectors) affords simple referencing in equations. 

As a second convention to minimize notation, for n > 
3, let the first row and first column of U be free of nega- 
tive signs in a CK parameterization, and let the top row 
be free of explicit complex conjugate symbols. 

For simplicity, let U G SU(n) be special unitary in 
this section, since given a special unitary matrix we need 
only multiply by a unit complex number to get a general 
unitary matrix, and other sub-groups can be reached via 
any additional general unitary transformations. 

Letting n be the dimension, there are two special cases. 
If n = 1, the only special unitary matrix is the scalar, 

U [1] = 1, (3) 

where numbers in square brackets of superscripts denote 
dimension n. For n = 2, U is not unique, and we have 

uin - ( 4- «• ) • < 4 > 

where a and b are complex numbers where \a\ 2 + \b\ 2 = 1. 



2 



General Formula for n > 1 



B. Examples for Cayley-Klein Form 



For all dimensions, n > 1, is given by 

UM = $ M -Q g^y-a;^^), (5) 

where products are written left to right, parenthetical 
subscripts on matrices indicate labels, not matrix ele- 
ments, and where fi^jj ^ (x, y) is the qubit factor matrix, 

(6) 

where 5j t k is the Kronecker delta and m is the sub- 
space identity matrix, 



(«,/3) 



y E\ n] 

/ j (m,r 



(a,/3) / ^ (m,m)> 

m— 1 



(7) 



in which ^j™'^ = \a)(b\ is the n-dimensional elementary 
matrix with a 1 in the row-a, column-b entry, and else- 
where. The matrix M^j^(s,t) from (6) is given by 



[n] , , = Ql,l(«, *) E (p, P ) + Ql,2(«, 4)^,0) 
JW (a,/3)V s > r J - n , Jn] n / , w n] 



(8) 

where a > /3, and Q(ui,«2) is the single qubit matrix, 

which has its alphabet of ordered variables reversed from 
that of (4) . The factor of 7„ in (6) is defined by 

Tn _E (_l)i(2— l + (-D"). (1Q) 

The matrix from (5) is the modified flip-matrix, 

n-l 

$h = 7))£ n (ii) 

/ / (m,n—m+l) V / 

m=l 

If 7„ = 1 in (11), $["1(7„ = 1) would oscillate between 
special and anti-special, such that starting at n = 1, 
det($M( 7n = 1)) = {1,-1,-1,1,1,-1,-1,...}, which 
is precisely the sequence generated by 7„. Thus j n en- 
sures that is always special without introducing any 
complex factors. The complex number pairs a a< p and 
b a ji from (5) are CK pairs, each subject to the constraint 

that |a Qj( g| 2 + |fc Q ,/3| 2 = 1. There are " ~" such pairs 
of parameters, with each pair corresponding to a spe- 
cial unitary transformation on a particular single qubit 
subspace, one corresponding to each lower-off-diagonal 
element, which are used to index the qubit rotations. 



To demonstrate the use of (5), for n = 3, if we let 

(02,1,62,1) = (e,/) 

(03,1,^3,1) = (a, 6) (03,2,63,2) = (c,d), 
then (5) produces 



(12) 



(a be bd 

b*e -a*ce-d*f* -a*de + c*f* |, (13) 
b*f -a*cf + d*e* -a*df-c*e* 

subject to the constraints \a\ 2 + \b\ 2 = 1, |c| 2 + \d\ 2 = 1, 
and |e| 2 + |/| 2 = 1. Notice that the first row and column 
conform to our conventions, and that in abbreviating the 
parameters in (12), we started with those corresponding 
to the bottom left corner, and then moved right along the 
bottom row, and then moved up the first column. For 
n = 4, if we define the abbreviations 

(02,1,62,1) = (j,fc) 

(03,1, 63,1) = {g, h) (a 3)2 , 63,2) = (I, m) 
(04,1,64,1) = (a, 6) (a 4)2 , 64,2) = (c,d) (a 4>3 , 64,3) = (e, /), 

(14) 

then, from (5) we obtain 



UW = 
( « 

6*3 



V 



be 

-a*cg 
+d*hl 



-a*ch*j 
-d*g*jl 
,+d*k*m* 



-a*ch*k 
-d*g*kl 
-d*j*m* 



subject to \a\ 2 + \b\ 2 

I2^|/j|2 = lf |j-|2 _j 



bde 
—a*deg 
—c*ehl 
+f*hm 
/-a*deh*'j\ 
'+c*eg*jl 1 
—c*ek*m' 

-f*g*jm 

\-f*k*l* 
f-a*deh*k\ 
+c*eg*kl 
+c*ej*m* 
-f*g*km 

\+f*j*r . 

+ \d\ 2 



bdf 
a*dfg 
c*fhl 
e*hm 
/-a*dfh*j\ 
+c*fg*jl 1 

-c* fk*rri 
+e*g*jm 
\+e*k*l* 
/-a*dfh*k\ 
1 +c*fg*kl ' 
+c* fj*m 
+e*g*km 
-e*j*l* 



1, 



(15) 

2 + l/| 2 = l, 



1. Sym- 



= 1, |c 

sr + i/ir = i, ur + |fc| 2 = 1, and |/| 

bolic and numerical tests up to n = 7 verify that (5) is 
special unitary and conforms to all proposed conventions. 
To reposition the single-term vectors of U^ n \ use 

U (r,c\r ,c )- ael l J (r ,r) D ( Co ,c)^(ro,r) L '(ro,co|ro,co)' ; '(co,c)' 

where (r, c) are the row and column to which you want 
the single-term vectors to move from (r , c ), and 



S n] ,= (1- 6 Z1 Z2 )E\ n] ,+ E n] ,+ I n] > 

(Zl,Z 2 ) ^ Z 1 ,Z 2 J (z U Z 2 ) (Z2,Zl) {Z1,Z 2 ) 



(17) 



is a swap matrix. The scalar in (16) keeps it special but 
adds an explicit phase factor which can be removed from 
the single-term vectors by redefining their Cayley-Klein 
pairs as needed. Note that (5) produces t/M = l/'"'^ ^. 



3 



III. HYPERSPHERICAL PARAMETERIZATION 
OF SPECIAL UNITARY MATRICES 

The Cayley-Klein parameterization convention devel- 
oped in the previous section gives us a methodical and 
orderly framework for describing special unitary matri- 
ces. Here, we adopt additional conventions which ensure 
a clear correspondence between the leading state vectors 
of U and a hyperspherical representation of a pure state. 



A. Review of Hyperspherical Coordinates 

An (n — l)-sphere is the surface of all points equidis- 
tant from a single point in n-dimcnsional Euclidean space 
(n-space). Therefore, in 3-space the 2-sphere is an ordi- 
nary sphere. An (n — l)-sphere can be parameterized 
by a single radial coordinate r, and n — 1 + 5 n _\ angular 
coordinates, n — 2 + 5 n> i of which are polar coordinates 
each with range [0, tt] , and one of which is azimuthal with 
range [0, 2it) if n > 2, (or with two-point range {0,7r} if 
n = 1), as shown for n > 2 below [9], 



Radial: 



Polar: 



Azimuthal: 



r e [0, oo) 
f 6i G [0,tt] 

?n-2 e [0, Tl] 

?„_i S [0,2tt) 



(18) 



For n > 1, a formula to express the n Cartesian coor- 
dinates of n-space using hyperspherical coordinates is 



/fc-i 

/ fc-2 



Xl<k<n- 
k=n^l 



1+6, 



Ar,{0u}) = r 



n 



Sill 







cos(0 fc ) 
sin(0 fe _ 



^m— 1 

(19) 

where the top row is the k th Cartesian coordinate where 
k 7^ n unless n = 1, and the bottom row is the n th 
Cartesian coordinate when n^l, empty products are 1, 
and {9 U } = {0i, . . . , 9 n -i + s n Inverses are 



r — \ / S i 

y m=l 

cos(0 n _i) 



9\<k<n-2 = cot 1 
■■, sin(0„_i) = 




valid for n > 3 for simplicity, though a more general 
form exists, and the azimuthal angle is left in its sine 
and cosine forms since both arc needed to find 9 n -i- 



B. Representing Pure States on Hyperspheres 

A pure state with global phase has the form 

(21) 



k=l 



with real moduli ri<k<n € [0,1], real phase angles 
4>\<k<n € [0,27r), and kets |1 < k < n) forming a com- 
plete basis in a Hilbert space of n dimensions, and 



k=i 



= i. 



(22) 



which suggests that the (non-negative) rfc can be de- 
scribed in the first sector of a unit hypersphere in n- 
space, while the <pk can be represented on n separate 
circles. Thus, for n > 1, the {/•&} = {n, . . . , r n } can be 
parameterized as 



r k = x k (l,{9 u }); u e(l-(y„,i)[O,f 



2J' 



(23) 



where Xfc(r, {^ M }) is given in (19), and each phase angle 
4>k is already its own parameterization. The mapping of 
the rfc in (23) only uses the first sector of a unit hyper- 
sphere because remapping to a full hypersphere does not 
give unique locations to all basis states. This method is 
in contrast to the Bloch sphere which uses relative phase 
together with the polar angle to put all coordinates into a 
single object. Here, the geometrical interpretation is that 
the superposition properties of \ip) are represented on a 
hypersphere, while the absolute phase of each outcome 
is represented on a separate phase circle, independent of 
the superposition hypersphere. 

Inverse equations for pure states when n > 1 are 



01<fc<ra-l+5 n ,i — cot 
4>i<k<n = arg(c fc ), 




(24) 



where c k = (k\ip) = r k e l ^ k , and arg(re^) = 4> e [0,2tt). 
An example of (23) applied to (21) for n = 4 is 



cos(6»i)e^ 1 |l) 
+ sin(0i)cos(6» 2 )e^ 2 |2) 
+ sin(0 1 )sin(0 2 )cos(0 3 )e l<fe |3) 
+ sin(0 1 )sin(6»2)sin(6'3)e^ 4 |4), 



(25) 



where 0i, 2 , #3 € [0, f] and 0i, 4> 2 , </> 3 , <j)A € [0, 2ir). Thus, 
the superposition properties of this state are portrayed 
by three angles 9 k in the first sector of a unit hypersphere 
in 4-space, while there are four separate absolute phase 
circles portrayed by the (j) k , or if global phase can be 
discarded, just three relative phase circles portrayed by 



»2|1 



»3|1 



and 



'4|1 



In general, a pure state has one superposition hyper- 
sphere and n — 1 relative phase circles, resulting in n 
separate geometric objects to describe a pure state in 
n-dimensional Hilbert space. An alternative geometrical 
approach is to use 2n — 1 Cartesian coordinates for the 
complex parts of (21) without global phase to represent a 
pure state on a single hypersphere in (2n— l)-space. How- 
ever, that method obscures the physical roles of its pa- 
rameters with respect to superposition and phase, while 
the method of (23) yields clear relationships between its 
parameters and physical properties. 
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C. Application to Unitary Matrices 

At first glance, it would seem that the most general 
way to build the CK constraints directly into f/M is 

{a J , k ,b j , k ) = (c j , k e^-';8 j , k e ie ^); ? < < " _ i, ( 26 ) 

where c jyk = cos(i? Jifc ), s jife = sin(i? Jife ), and i?^, (/) jtk ,e jtk e 
[0, 27t). However, this introduces more variables than the 
true number D of degrees of freedom in It turns out 
that there exists an essential parameterization such that 
D for special and general unitary matrices is 

D S v(n) =n 2 -I and D u(n) = n 2 . (27) 

Since the parameterization of (26) uses 3( - 2 ~ n ) vari- 
ables, this means that it introduces — x ^ n ^ irrelevant 
parameters, a waste of computation time that scales as 
for large n. Therefore, to further constrain (26) to 
have only the essential parameters, let 

£j,k = Xj,k$j-i,k, (28) 

where \ G [0, 27r). Thus, (28) ensures that (26) only in- 
jects n 2 — 1 variables into U^. The meaning of this is 
that two phase angles are needed by all single qubit rota- 
tions whose nontrivial off-diagonal elements are directly 
next to the main diagonal elements, while all other ro- 
tations only need one phase angle. In terms of quan- 

2 

turn states, this introduces " 2 " superposition angles 
and " 2 +" — 1 phase angles for SU(n). To map the super- 
position angles onto first sectors of hyperspheres, use 

^, k = e jtk e(l-S nA )[0,^]. (29) 

Thus (29) allows us to relate the form of f/M to quantum 
pure states in a way particularly evident in its leading 
vectors. To obtain a general unitary matrix € U(n), 
simply apply the essential parameterization to , and 
then ict GN = e^UW, where 7 e [0,2tt). 

D. Examples of Hyperspherical Unitary Matrices 

To illustrate the benefits of the above parameteriza- 
tions, for n = 3, using (26) and (28) in (13) yields 

urn = 

/ T?4>3,l+T,0 p<^3,2,0 pX3,2,0 \ 

W 3,1,U,U 03,1" 2>»3,2,7T »3, 1-2" i»3,2- 2 

/ — 03,1+03,2 + 02,1,^ /-03,1+X3,2+02,1,\ 
P<l>2,lfl ^-\X3,2+X2,1 /p\03,2+X2,l / 

*3,l-f ,7T,1?2,1 03,1,03,2,02,1 03, 1 ,03,2 - f ,02,1 

/-03,l+03,2+X2,l,\ /^03,1+X3,2+X2,l"!\ 
^X2,l,0 ^Ax3,2+02,1 ^103,2+02,1 ^ 

V 03,1-f ,T,02,l-§ 03, 1,03, 2,02, l-§ 03,1,03,2-f >02,l-§/ 

(30) 

where, given c(x) = cos(x) and s(x) = sin(x), we define 
FS;* iX = ~c{v)c{w)c{xyy - s(w)s(x)e- iz . (31) 



Thus, a single function can be used to describe all ele- 
ments. Employing the first-sector hypersphere mapping 
of (29) then reveals the columnized first row of (30) as 

ri ( U i\\ / c(0 3 ,i)e^ \ 
PS) = \u\ 3 l = s{e^)c{^a)^ 3 ' 2 > (32) 
\U%) \s{h,M9 3 , 2 )jx™ J 

where 6 3 . 1 ,e 3 . 2 G [0, f], and 03,1, ^3,2, X3,2 € [0, 2?r), so 
this vector has the standard hyperspherical form of a pure 
state, with a similar result holding for the first column 
of , in keeping with our conventions. 

For n = 4, the first row of has the form of (25), 
and using the methods presented in this paper, this same 
pattern is maintained for all dimensions n > 1. The true 
benefit of the form in (30), and of this method for all n, 
is that it incorporates the CK constraints automatically, 
while using the fewest parameters possible. 



IV. PHYSICAL APPLICATIONS 

While the conventions presented in this work arc in- 
teresting mathematically, their worth to physicists will 
become apparent in the many ways their use can stream- 
line and enhance the field of quantum mechanics. To il- 
lustrate this, we now discuss a few physical applications. 



A. Complete Basis Generation 

Given a state vector \?p) in a Hilbert space of dimension 
n > 2, suppose we wish to find a set of n — 1 vectors that 
are orthonormal to it. The standard prescription for this 
task is the Gram-Schmidt orthonormalization procedure, 
but it can be tedious to use and it requires that we start 
with a set of n vectors. 

An interesting alternative arises from the methods of 
this paper. First, note that for a unitary matrix con- 
structed from (5), the leading vectors have single-term el- 
ements. This yields simple expressions for the first 2n — 1 
essential parameters of a unitary matrix by using the hy- 
perspherical inverse equations of in (24). Then, we 
use the property that all the rows of a unitary matrix 
form a complete basis, obtaining our n — 1 new vectors 
from the remaining rows of up to n(n~ 2) undeter- 
mined degrees of freedom that can be freely chosen. 

Therefore, to find the complete basis of states that 
includes the known vector, 

n 

foft) = (33) 
k=i 

first get the superposition and phase angles of (33) from 
(24). Then, to relate them to individual qubit rota- 
tions using our conventions, use the ordering described 
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in the Cayley-Klein examples to relabel these angles as 
the dual-indexed angles, 

{O n ,k,(j>n,k) ={9 k ,4) k ); \<k<n-l ^ 

Xn.n— 1 — 071- 

Thus, defining the CK pairs using (26), (28), and (29), 
all essential parameters whose first index is n have val- 
ues given by (34) which defines the determined CK pairs 

( a nA, ^n,l)i ( a n,2, 6/1,2), • • • , (o n ,„_i, 6n,n-l)- 

The remaining CK pairs can be freely chosen. Specifi- 
cally, we can arbitrarily choose all (oj,fc, bj, k ) for 2 < j < 
n — 1 and 1 < k < j — 1 such that |flSj,fc| 2 + |6j,/c| 2 = 1- 
These constraints can be satisfied automatically with the 
fewest parameters by using the essential parameteriza- 
tion of (26), (28), and (29), which then requires only 
n(n — 2) free parameters to be chosen. 

Then, putting the determined and freely chosen CK 
pairs into (5) produces a matrix t/M where the first row 
vector is the generating state and the remaining n — 
1 row vectors are orthonormal to \ip) and each other. 
Therefore, using the method just described, it is possible 
to generate a complete basis of states from a single input 
state \ip) and n(n — 2) free real scalars. 

For example, to find an orthonormal basis that includes 

\1>) = ci|l> +c 2 |2) +c 3 |3) +c 4 |4>, (35) 
first use (24) to get its hypcrsphcrical coordinates as 

^-^i^'h-^M- (36) 
9 3 = cot and 0i< fc < 4 = arg(ci< fe < 4 ), 

where r k = \ck\ = |(fc|V>)|- Then using (36) in (34) gives 

(04,l,04,l) = (01, 0l), (#4,2,04,2) = (02,02), (04,3,04,3) = 

(03,03), and X4,3 = 4 , which then determine the essen- 
tially parameterized CK pairs (014,1, 64,1), (04,2, 64,2), and 
(04,3,64,3). Using the abbreviations of (14), this gives us 
the determined CK pairs as 

(0,6) = (04,1,64.1) = (cos(0i)e*,sin(0i)) 

(c,d) = (04,2,64.2) = (cos(0 2 )e^,sin(0 2 )) (37) 

(e,/) == (04,3,64,3) = (cos(0 3 )e^,sin(0 3 )e^), 

where these 7 angles are found in (36). The free CK pairs 
are then given in essential parameterization as 

(g,h) ee (03.1,^3.1) = (cos(^3,i)e^ 3 - 1 ,sin(^ 3 ,i)) 
(j,k) = (o 2 .i,6 2 .i) = (cos(^ 2 ,i)e^. 1 ,sin(fe,i)e^. 1 ) 
(Z,m) = (03.2,63.2) = (cos(^3,2)e^ 3 . 2 ,sin(^, 2 )e^3,2) ; 

(38) 

where the 8 free parameters 02,1,03,1,03,2 € [0, § ] and 
02,1, 03,1, 03,2, X2,i,X3,2 € [0, 2tt) can have any values. 
Then, using (37) and (38) to generate 17 M from (5) as 
shown in (15), using notation similar to (32) we obtain a 
complete set of orthonormal state vectors, given by 

{mu%b,\41.),\uS)}, 09) 

which are the transposed row vectors of where the 
generating state is \ip) = \U^}). 



B. Multiport Interferometry 

As described in [8], a multiport interferometer (MPI) 
can be constructed by performing successive single qubit 
rotations on all possible pairs of spatially separated 
beams. The resulting operation is a general unitary ma- 
trix B of dimension n, the number of input ports. 

In such a setup, one uses a two-arm Mach-Zhender in- 
terferometer to implement each single qubit operation, 
with phase-shifter angles acting as the superposition and 
phase angles of each rotation. As described earlier, if 
each qubit rotation is general, there are a large number 
of irrelevant parameters. Therefore, using the essential 
parameterization method of (26) and (28), one can elim- 
inate unnecessary hardware to achieve the simplest ex- 
perimental setup. Furthermore, by also using the first 
sector mapping of (29), one can minimize the angular 
ranges required from the superposition phase shifters. 

Another benefit of these methods stems from the fact 
that the output spatial field modes of the MPI are 

n n 

a' k ee Ba k B^ = ^ V k , m a m = ^2 U m,k^m, (40) 

m— 1 m— 1 

where a m is the annihilation operator for the input field 
mode of location m, primes indicate output modes, and 
V and U are general unitary where U = e^U^ = V\ 
and U will simplify what follows. Note that V and U arc 
not the same as the unitary B of the MPI, which can be 
viewed as a multiport variable beam-splitter. However, 
the transformation of (40) can be expressed directly in 
terms of phase shifter angles, and can be derived for a 
particular set of equipment using the methods of [8] . 

An example of a particular application is the explicit 
calculation of the output of an MPI given the input of a 
product of different coherent states. By using (40) and 
the fact that coherent states are \a) = D(a)\0) where 
D(a) ee e Qat - Q * a i s the displacement operator, we find 

B ® \a k ) = | ELi U k , m a m ) ee ® \a' k ), (41) 

k=l k=l k=l 

which yields the output as a product of coherent states, 
and also allows us to use the compact CK form of (5) to 
describe the new coherent state parameters (CSPs) a' k . 

For example, given a 3-input interferometer, with three 
coherent states \a), |/3), and I7) in product as input, (40) 
and (41) reveal the MPI output to be 

B\a)\f3)\ 7 ) = \a')\P')\ 7 '), (42) 

where, using U = e^U^ from (13), the new CSPs are 

a' ee e^(aa + bc/3 + bdry) 

/3' ee e^(b*ea + (-a*ce - d*f*)/3 + {-a*de + c* f*)j) 
i ee e**(6*/a + {-a*cf + d*e*)/3 + {-a*df - c*e*) 1 ). 

(43) 

Thus, plugging the essential parameterization of (28) into 
(43), with the help of (12), allows us to determine pre- 
cisely how the phase shifters control the coherent state 
parameters in the output state of (42). 
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From (43), the sum of square magnitudes of the CSPs 
is invariant; \a'\ 2 + |/3'| 2 + | 7 '| 2 = \a\ 2 + |/?| 2 + | 7 | 2 , so 

which is because the vectors of U arc orthonormal. 

A physical application of all this is the ability to cal- 
ibrate a multiport interferometer with classical inputs. 
By setting different combinations of the coherent state 
inputs to vacuum in (42), one can isolate any problems 
in the output of a large MPI to only a few possible cul- 
prit qubit rotations. Then, once the proper U is verified 
using the convenient classical calibration beams, one can 
remove them and proceed with any other input desired. 

C. Entanglement Preserving Transformations 

The methods of this paper are also useful in finding ex- 
plicit forms for entanglement preserving transformations. 
For TV-party systems, these are typically generated as 

■ «1 \[nil . CIN \\n*r] 

U EP = e r ' , e- , ~2- u < 1 >- A ® ■ ■ ■ ® e -^ u <«>- A ' Nl , (45) 

where 7 and a m are angles, the U( m ) are real unit vectors 
of dimension n m , and X^ Um ^ are vectors of generalized 
Gell-Mann matrices of dimension n m . However, matrix 
exponentiation can be avoided completely if we use (5) 
with the essential parameterization of (28) to write 

f/fip-e 17 ^ 11 ®-®^; 1 , (46) 

where subscripts k m indicate that each unitary matrix 
has a separate set of essential parameters, resulting in 
Dep = 1 — N + X^m=i n m variables, the same as in (45). 
For example, in a qubit-qutrit system, if 7 = in (46), 

U E p = Uf ] ® f/f = 

fa b\ ( 0,2 b2 ° 2 &2 ^ 2 \ 

\-b* a* ) ® &262 _a 2 C2e 2 - ^2/2 -a* 2 d 2 e 2 + c^/ 2 * , 
V 1 lJ V&2/2 -a* 2 c 2 f 2 + d*e* -a*d 2 f 2 - c*e* 2 J 

(47) 



which is subject to (26) and (28). Thus, the essential 
parameterization lets us express entanglement preserving 
transformations without the need for exponentiation. 



V. CONCLUSIONS 

The methods presented here produce unitary matrices 
in an orderly, predictable fashion, highlighting the con- 
nection between unitary matrices, pure states, and hy- 
perspherical coordinates. These conventions reveal the 
true structure of unitary matrices with constraints al- 
ready built-in, simplifying tasks such as generation of or- 
thonormal bases, multiport intcrfcrometry, and generat- 
ing entanglement preserving transformations. The meth- 
ods of this paper may also be useful for differentiation, 
optimization, approximation, and quantum control. 

The presented conventions are designed to be conve- 
nient for many applications in quantum mechanics, but 
they are not unique. Over the years, several methods 
have been developed to parameterize unitary matrices, 
such as in [10-15], and some methods may be more use- 
ful than others depending on the application. 

The most general unitary matrix, as expressed using 
the conventions of this paper, is G^ = e n U^ n \ where 
7 G [0, 27r), and £/N is given by (5) and parameterized 
by (26) subjected to (28) to ensure that Gl"l has ex- 
actly n 2 variables. Thus, this method allows us to con- 
struct the most general unitary matrices using the fewest 
parameters possible without the need for exponentiating 
combinations of generators of U(n), producing a simple 
form that is both predictable, and readily applicable to 
a wide variety of fields of interest in physics. 
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